Dynamics of fluctuations below a stationary bifurcation to electroconvection in the 

planar nematic liquid crystal N4 
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We fitted C(k, t, e) oc exp[— <r(k, e)r] to time-correlation functions C(k, r, e) of structure fac- 
tors S(k, t, e) of shadowgraph images of fluctuations below a supercritical bifurcation at Vo = 
V c to electro-convection of a planar nematic liquid crystal in the presence of a voltage V = 
\/2Vocos(27r/t)[k = (p, g) is the wave vector and e = V 2 /K, 2 — 1]. There were stationary oblique (nor- 
mal) rolls at small (large) /. Fits of a modified Swift-Hohenberg form to <r(k, e) gave /-dependent 
critical behavior for the minimum decay rates ao(e) and the correlations lengths £, P , q (e)- 

PACS numbers: 05.70.Jk, 05.40.-a, 64.60.Fr, 45.70.Qj 



Critical opalescence is a familiar phenomenon that oc- 
curs, for instance, near liquid-gas critical points and near 
the phase separation of binary mixtures. [l| It is caused 
by refractive-index fluctuations due to droplets of fluid 
of lesser or greater than average density or concentra- 
tion. The fluctuations are the response of the system 
to the thermal noise inherent in the Brownian motion 
of the molecules or atoms. The droplet size increases in 
proportion to a correlation length £ that diverges at the 
critical point. Likewise, the decay rate a of the droplets 
vanishes as the critical point is approached. When £ 
is sufficiently large, ambient light is scattered and the 
fluid becomes "milky" . Analogous phenomena occur in 
many spatially extended nonlinear systems as they are 
driven away from equilibrium and toward an instability, 
or bifurcation. |3,y,y,|5|, |a, |7|,|3| For those cases the length 
and time scales are more nearly macroscopic and the fluc- 
tuations can be observed directly and in detail with mod- 
ern shadowgraph techniques. 0, [ToL 1 1 1 1] Here we present 
experimental results for £ and a for electro-convection 
(EC) of a nematic liquid crystal (NLC). As we will show 
below, our data are not easily understood by analogy to 
known equilibrium critical phenomena Il2| . 

A NLC consists of elongated molecules that tend to 
align locally relative to each other. The alignment 
direction is called the director n. When a NLC is con- 
fined between parallel glass plates with a small spac- 
ing d between them, n is influenced by the interaction 
of the molecules with the glass surfaces. We used sur- 
faces that caused more or less uniform alignment in a 
unique direction parallel to the plates (planar alignment). 
This breaks the rotational invariance usually found in 
isotropi c liq uids, for instance in Rayleigh-Benard convec- 
tion, pf [lj, |l5j or in NLC with homeotropic alignment 
(n perpendicular to the plates) [l^.We applied a volt- 
age V = V2V~o cos(2n ft) between transparent indium-tin 
oxide electrodes on the inner surfaces of the confining 
plates. For a NLC with suitable properties this induces 
a transition (or bifurcation) at Vo — V c from a spatially 
uniform state without flow to a convecting state of lower 



symmetry. |17| The bifurcation is supercritical, analo- 
gous to a critical point in a two-dimensional equilibrium 
system, e.g. to a Curie point of a ferromagnet. Here we 
report on the effect of fluctuations on the nature of the 
bifurcation to EC. 

When the fluctuations are small, one may neglect in- 
teractions between them and discuss t he p henomena in 
terms of a linear theory (LT).[l3. HsL Hfl| In that case 
we expect that £ ~ N" with v = 1/2 and a ~ |e| A 
with A = 1. Here e = V<? /V? - 1. When the fluctu- 
ations become larger, nonlinear terms in the governing 
equations lead to interactions between them. In that 
case several universality classes of critical behavior ex- 
ist for equilibrium systems. They are determined pri- 
marily by the dimensionality of the system and by the 
number of components n of the order parameter. For 
two-dimensional equilibrium systems one has primarily 
Ising systems (n = 1) and X — Y systems (n = 2). It is 
not certain that this classification carries over in direct 
analogy to non-equilibrium cases. 

The range of e over which fluctuation interactions be- 
come important depends on the coupling of the thermal 
noise to the system. For EC, the relevant parameter is 
the ratio of the thermal energy to a characteristic elastic 
energy of deformation F = ksT/kd (k is an average elas- 
tic constant of the NLC).Q On the basis of a Ginzburg- 
like criterion!^ one can estimate that the critical region 
is encountered near e c ~ F 2 / 3 . For the NLC used by 
us one finds k ~ 10~ n N [2l| and thus, for a cell with 
d ~ 25 fiva, one has e c ~6x 10~ 4 . For much larger e the 
prediction of linear theory should pertain. 

We made shadowgraph measurements at various values 
of /. We note that a priori there is nothing to indicate 
that e c or the universality class should depend on /. We 
used the same 0.88 x 0.88 mm 2 area of the same sample 
for all measurements, thus reducing the possibility that a 
frequency dependent result could actually be due to dif- 
ferent sample preparation. As shown in Fig.^i and c, at 
large / = 4 kHz we found so-called "normal" stationary 
convection rolls in the fluctuations below and the pat- 
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terns above onset. Their wave vectors are parallel to n. 
For the fluctuations we found A ~ 1, consistent with LT. 
However, the correlation lengths gave v — 0.40 ± 0.03, 
which differs from v = 1/2 of LT. One might argue that 
this measured exponent really is an effective exponent 
representing data in a crossover region from linear to 
critical behavior. However, then one would have to ar- 
gue either that the crossover region is different for £ and 
a or that the value of A is not altered measurably by 
fluctuation interactions. 

At small / = 25 Hz we find that the fluctuations, 
as well as the patterns above onset, correspond to 
oblique stationary rolls with their wave vectors forming 
an oblique angle \6\ > with n (Fig.^ 3 an d d). For this 
case there are two degenerate modes, known as "zig" 
and "zag" , with angles ±0. Here the critical behavior 
was qualitatively different. No power law could describe 
the results for cr and £. Both quantities, although they 
initially changed upon approaching V c , deviated from a 
power law for |e| < 0.02 and approached a constant value 
rather than diverging. We find it difficult to reconcile 
these results with any known critical behavior of equilib- 
rium systems. 

We took shadowgraph images 0, 0, 0] I, (x, e) [x = 
(x, y) are the spatial coordinates] of a commercial cell [22] 
with a nominal spacing of d = 25 and filled with the 
NLC Merck N4 (a eutectic mixture of two azoxy com- 
pounds, CH30-C6H4-NO=N-C6H4-C4H9 and CH30- 
C6H4-N=NO-C6H4-C4H9) doped with 0.1% by weight of 
tetra butylammonium bromide. The alignment was pla- 
nar. The conductance was measured at 30°C, at / = 40 
Hz and V = 2.0 Volt. It was near 6.1 x 10 _7 (flm) -1 . We 
used the total shadowgraph power P to determine that 
the bifurcation to EC was non-hysteretic and stationary 
over the range 25 < / < 8000 Hz. The power above V c 
was used to find V c = 6.04 (11.49) Volt for / = 25 (4000) 
Hz. Above the bifurcation the pattern consisted of sta- 
tionary zig and zag (normal) rolls for / < /x, ~ 1.3 kHz 
(/ > /l)- At each voltage we waited 300 s, and then took 
4096 shadowgraph images 0.3 s to 1 s apart, depending 
on e. 

From each image we calculated Ii (x, e) = 
Jj (x, e)/I(j (x, e) — 1. Here 2o(x, e) is a background 
image obtained by averaging 4096 images at the same 
e. For each ij(x, e) we derived the structure factor (the 
square of the modulus of the Fourier transform) Si (k, e) 
and averaged 4096 $j(k, e) to get averaged structure 
factors S*(k, e), where k = (p, q) is the wave vector. 
Figure ^c) and (d) show examples. At low frequencies 
[(d), / = 25 Hz] two pairs of peaks correspond to 
two sets of rolls oriented obliquely to the director. At 
high frequency [(c), / = 4 kHz] only one pair of peaks 
corresponding to normal rolls was observed. 

Figure [3 shows time series of 5(ko, e, t) for e = —0.010 
at k = ko where S^k, e) has a maximum. The upper 
(lower) figure shows the result at low (high) frequency. 
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FIG. 1: (a,b): Shadowgraph images of size 0.88 x 0.88 mm 2 
for e = -1.1 xlO" 3 . (a): / = 40QQHz. (b): / = 25Hz. (c, d): 
the central part — 7.2 < p,q < 7.2 of the averaged structure 
factor S(k) based on 4096 images for e = -1.1 x 10" 3 . (c): 
/ = 4000#z. (d): / = 25Hz. The director is horizontal. 
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FIG. 2: Time series of S(kj>,(?)) at e = -1.0 x 10" 2 for 
(a) / = 25 Hz, k = (3.936, 1.968) and (b) / = 4000 Hz, 
k = (4.830,0). 

The fluctuations of the structure factor are random and 
of similar size at both frequencies. 

In Fig. we show the auto-correlation function 
C(k,r) =< S(M)S(k,t + t) > I < S(k,t)S(k,t) > 
(< ... > indicates an average over t) of the fluctuations 
shown in Fig. 0[a) on linear scales. The insert gives 
the same data on a semi- logarithmic plot. One sees that 
the correlations decay exponentially in time. A fit of 
Ae~ CTT to the data yields a = 0.14 s^ 1 . Although S(k) 
differs from the structure factor of the refractive index 
by a fc-dependent factor [the shadowgraph transfer func- 
tion T(/c)],[ll| we note that C(k, r), and thus a, are not 
influenced by T(fc). 

Decay rates of fluctuations at £ = —0.010 are shown 
in Fig. 0] for oblique (a, c, / = 25 Hz) and normal (b, 
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FIG. 3: The time auto-correlation function of the fluctuations 
of the structure factor at k = (3.936, 1.968) for / = 25 Hz 
and e = -1.0 x 10" 2 . 



d, / = 4 kHz) rolls along two lines parallel to the p and 
q axes and passing through the maxima of S(k) at ko. 
In the p (q) direction and for / = 25 Hz a is shown in 
Fig. 0k (c). The Swift-Hohenberg forms 



a(p,q ) = cro(e) x [(e p ( P 2 - p 2 ) 2 + 1} 



with qo — 1.968 and 



<r(Po,q) = <ro(e) x [£(q 2 



(1) 



(2) 



withpo = 3.757 were fitted to the data (solid lines). From 
Fig. QJ: it is clear that a simple quadratic form, restricted 
to either positive or negative q, would not provide a good 
fit. 

The large-/ results for a are shown in Figs.0J) and d 
for the p and q direction respectively. Again, we fitted 
Eq. [I] to the data in Fig. (solid line). However, for 
normal rolls we fitted the function I23J 



a(p = 0,q)= a (e) x (^g 2 + t\d + 1) 



(3) 



to the data (solid line). 

Each fit yielded a value at a given e of the minimum 
decay rate ao(e), an d of the correlation length £ p (e) or 
£q(e). For / = 4 kHz we plot £ p and cr as a function 
of |e| in Fig. [3^ and b respectively. Fits of the results 
for |e| < 0.05 to power laws yielded v = 0.40 ± 0.03 and 
A = 0.92 ±0.06. The result for A is reasonably consistent 
with the prediction of linear theory and suggests either 
that the critical region has not been reached in this ex- 
periment, or that fluctuation interactions do not alter A 
by a measurable amount. The result for v differs sig- 
nificantly from the LT prediction. Of course one could 
argue in this case that the values of |e| are too large for 
LT to apply; but in that case the range of applicability 
would have to be different for £ p and a . Thus the data 
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FIG. 4: The decay rates o(p,q,e) as a function of p and q 
at e = -1.0 x 10~ 2 for (a,c) / = 25 Hz and (b,d) / = 4000 
Hz. The solid lines are fits to the data that yielded (a) ao = 
0.13, = 0.37, (b) a = 0.21, £ p = 0.29, (c) a = 0.18,6, = 
0.21, and (d) a = 0.31, £ g2 = 0.87, £, 4 = 0.11. 
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FIG. 5: (a) The correlation length £ p and (b) the decay rate 
cto versus |e| for / = 4 kHz derived from fits of Eq. to 
<y(p,q = 0). The dashed line in (a) corresponds to v — 1/2. 
The solid lines are fits to the data that gave (a) v = 0.40±0.03 
and (b) A = 0.92 ± 0.06. 



suggest, but do not firmly establish, a deviation of the 
critical behavior from the LT prediction. 

In Figs. EJi and b we show analogous results for / = 
25 Hz. Here the situation is quite different. For |e| < 
0.02 deviations from the LT prediction are noticeable for 
do, £p, and £ q . Indeed, a power law can not be used to fit 
the data. The results suggest that all three parameters 
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FIG. 6: (a) The correlation lengths £ p (solid circles) and f 9 
(open circles), and (b) the decay rates oo derived from fits 
of Eq.Qto u(p,qo) (open circles) and Eq.|5]to a(po,q) (solid 
circles), versus |e| for / = 25 Hz. The solid lines correspond 
to v = 1/2 and A = 1 in (a) and (b) respectively.. 

approach a finite value at e = 0. This behavior differs 
from that of any known equilibrium system near a critical 
point. However, it would be consistent with a first-order 
phase transition in the presence of strong noise. 

We also measured the total power P = / n °° S(k)cPk at 
various |e|, and found that the |e| dependence of P de- 
pended on /. Correlation lengths derived from the widths 
of S(k) agreed with those reported here. At small / we 
also found that the maximum of S(k) remained finite, 
whereas one would have expected it to diverge at e = 0. 
Those results also differed from previous measurements 
for EC near a Hopf bifurcation in the NLC 152. [101 This 
work will be reported in a more detailed publication. 

In this Letter we presented measurements of the dy- 
namics of thermally driven fluctuations near the bifurca- 
tion to stationary electro-convection of the nematic liq- 
uid crystal N4. The results for the decay rate Co and the 
correlation lengths £ p and £ g for oblique rolls at small 
driving frequency / tend toward a constant value at the 
bifurcation. This differs from the prediction of linear 
theory which yields power-law singularities with expo- 
nents A = I and v = 1/2 respectively. It also differs 
from known possibly relevant equilibrium critical behav- 
ior which yields vanishing decay rates and diverging cor- 
relation lengths, albeit with modified exponent values. It 
is consistent with a first-order transition in the presence 
of strong noise where no hysteresis would be observed. 
The results for normal rolls at large / are more nearly 
consistent with LT, although the value of v is somewhat 
lower than the predicted value 1/2. There was no a pri- 
ori reason to expect a dependence of the results on the 



driving frequency. 
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